Prom multiplicative noise to directed percolation in wetting transitions 
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A simple one-dimensional microscopic model of the depinning transition of an interface from an 
attractive hard wall is introduced and investigated. Upon varying a control parameter, the critical 
behaviour observed along the transition line changes from a directed-percolation to a multiplicative- 
noise type. Numerical simulations allow for a quantitative study of the multicritical point separating 
the two regions, Mean-field arguments and the mapping on a yet simpler model provide some further 
insight on the overall scenario. 
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A variety of interesting physical phenomena corre- 
sponds to the unbinding transition of an interface from 
a fiat surface. This is the case of wetting processes 
(WP) taking place in the thin liquid film which forms 
on a substrate exposed to a gas. By varying external 
parameters such as temperature or pressure, the liquid 
layer h(x, t) may either grow and become macroscopi- 
cally thick or remain confined to the close vicinity of 
the substrate 0, Q • Wetting phenomena can also take 
place under non-equilibrium conditions. Here one is in- 
terested, for example, in a growth process of a film over 
a substrate. Depending on the dynamical rates control- 
ling the growth process one can observe similar pinned 
or unpinned phases. 

A question of general interest concerns the universality 
of the unbinding transition. While the equilibrium sce- 
nario is well established fjj, an overall understanding of 
non-equilibrium wetting phenomena is still lacking. Nu- 
merical studies of non-equilibrium systems have revealed 
a composite picture that still needs to be fully disen- 
tangled. Here one considers the dynamical equations of 
a moving interface interacting with a hard wall. In the 
simplest case, the unbinding transition is signalled by the 
change of sign of the average velocity of the free inter- 
face. A pinned phase is obtained when the free interface 
moves towards the substrate and an unpinned phase is 
found when it moves away. This scenario is sometimes 
referred to as the multiplicative noise (MN) and is well 
described by a Kardar-Parisi-Zhang (KPZ) equation Q 
with a hard- wall 0, 0, IS • On the other hand, a study 
of some discrete growth models has shown that for partic- 
ular values of the growth rates the unbinding transition is 
of different nature, belonging to the directed percolation 
universality class |g. 

More interesting is the scenario when the surface ex- 
erts, in addition, a short range attractive force that may 
prevent h from growing even when the "free" interface 
would otherwise have a positive velocity. In WP such an 
effective interaction is obtained when the growth rate on 
the bare substrate is lower than the growth rate on the 



film itself. Some microscopic models, introduced to study 
this out-of-equilibrium depinning transition 0, , in- 
dicate that for a sufficiently strong attractive force, the 
unbinding transition may become discontinuous. 

This scenario is by no means restricted to WP. Similar 
features have been, e.g., observed in the onset of c omp lete 
synchronization in chains of coupled oscillators|lll IT^ . 
Il3j . In this context, two replicas of the chain are either 
subject to the same stochastic forcing, or locally coupled: 
the absolute difference A(i, t) between the state variables 
in the ith site may either go exponentially to zero, in 
which case synchronization eventually occurs, or stay fi- 
nite. Upon formally introducing j(i,t) = — logA(i, t), 
it is easily recognized that f) plays the same role 
as h(x,t): synchronization corresponds to the unpinned 
phase, while the unsynchronized regime corresponds to 
the bounded phase. In this context, a strong nonlin- 
earity may play the role of the attractive force in WP, 
preventing small but finite perturbations from vanishing. 
However, at variance with WP, in complete synchroniza- 
tion, the effect of a sufficiently strong "attractive force" 
is to bring the MN transition into the directed percola- 
tion (DP) universality class [Til IT^ rather than making 
it first order. 

In the absence of a sufficiently general field-theoretic 
approach able to reconcile all the various observed sce- 
narios into a common framework, the study of minimal 
models is very helpful for the identification of the basic 
mechanisms. This is the main motivation for introducing 
hereafter a simple microscopic model. By numerically re- 
constructing its phase diagram, we shall be able to clearly 
recognize that both the MN and DP universality classes 
can be found and to quantitatively investigate the "mul- 
ticritical" point separating the two scenarios. Moreover, 
an accurate reconstruction of the critical line will be pro- 
posed based part on mean-field arguments and exploiting 
the exact mapping (in a limit case) onto a genuine DP- 
model. 

Specifically, we consider a simple growth model, which, 
for reasons that will become clear in a while, is called 
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FIG. 1: Updating rule of the SSW model. The full line 
represents the interface, while the shaded area identifies the 
wall. Dashed segments indicate interface flips occurring in 
randomly chosen local minima (see A and B in panel a) and 
in all sites located below the wall after it has been shifted 
upwards by one unit (see C in panel b). 



Single-Step-plus-Wall (SSW) model. The starting point 
is the Single-Step (SS) model originally introduced to 
study the roughening of Id interfaces OH!!- It is well 
known that the SS model can be exactly mapped onto the 
Id KPZ equation^. The interface is described by a set 
of integer heights hi on the sites i of a one-dimensional 
lattice of length L, satisfying the "continuity" constraint 
\hi — h i+1 \ = 1. At each time step, dt = 1/L a site 
i is randomly selected and its height increased, hi — > 
hi + 2, provided that a local minimum exists at site i. In 
the thermodynamic limit, since the dynamics does not 
introduce any spatial correlation, a generic interface with 
mean slope s moves with a mean velocity v(s) — (1 — 
s 2 )/2. The exact knowledge of v(s) will be crucial in the 
following, since it allows for an exact determination of 
the critical line in the MN regime. 

The second ingredient of the SSW model is an upward- 
moving wall located at some integer height h w (t) below 
the SS interface. It moves with velocity v w and both 
"pushes" and attracts the interface. Altogether, the SSW 
dynamics amounts to the following evolution algorithm: 
at each time step, a site i is randomly chosen and, if it 
is a local minimum, hi is increased by two units with 
probability 1, or (1 — q), depending whether hi > h w , or 
hi = h w (see Fig. QJi). After n w — L/v w steps, the wall 
is moved upwards by one unit and, simultaneously, the 
height of the interfacial sites overtaken by the wall, is in- 
creased by two units (see Fig.^ 3 ). Physically, the SSW 
describes a roughening and moving interface, attracted 
by a short-range force to a hard wall. Its dynamics is de- 
termined by two parameters: (i) the relative velocity of 
the wall with respect to the free-interface, that we control 
by modifying 1^,0]; an d (ii) the stickiness of the wall, 
quantified by q. Since we are interested in characteriz- 
ing the phase diagram of SSW by locating the depinning 
transition from the wall, the natural order parameter is 



the density of sites pinned at the wall 

(i) 

where (•} denotes an ensemble average over different re- 
alizations of the stochastic process. 

A necessary condition for interface depinning to occur 
is v w < 1/2, because 1/2 is the velocity of a free and 
flat interface. For q small enough, this is also a sufficient 
condition, since the attractive force is overcompensated 
by the faster velocity of an almost straight interface. This 
can be seen by a simple mean field argument: an interface 
completely attached to the wall has a higher density of 
minima (1/2) than a rough one (1/4), so that its average 
velocity v = (1 — q) is larger than 1/2 as long as q < q* = 
1/2. Accordingly, below q* , a pinned interface detaches 
as soon as v w < 1/2. Numerical simulations confirm 
that the transition indeed occurs at v w — 1/2, with the 
only slight difference that q* — 0.4445(5). Above q* the 
interface may remain pinned even when its velocity is 
larger than v w (see Fig. 0). In other words there is a 
sort of bistable region, where an initially pinned interface 
remains attached while a depinned one moves away from 
the wall. The transition line, located at the lower border 
of the bistable region, is continuous and turns out to 
belong to the DP universality class. This is at variance 
with the discontinuous transition observed, e.g., in the 
solid-on-solid model of [T^j. Before commenting on the 
possible reason of such a difference, it is necessary to 
explore in more quantitative way the critical behavior 
both above and below q* . 

Continuous non-equilibrium phase transitions are char- 
acterized by three independent critical exponents. At 
criticality (v w = v^), the density pit) of pinned sites 
scales with time as p ~ t~ s , while its stationary value 
depends on the the distance from criticality as, 

limp(t)~K-<y, (2) 

t — >oo 

Very accurate numerical estimates of the DP critical ex- 
ponents give P DP = 0.276486 ± 6 • 10~ 6 and S DP = 
0.159464 ± 6 ■ 10 -6 ^?|]- L ess accurate estimates are 
available for the MN scenario, namely (3m n = 1.7 ± 0.1 
and Smn = 1-1 ± 0.1 5]. The third exponent, z, 
can be defined with reference to the scaling relation 
p(t,L) = L~ Sz g(tL~ z ) 0, where g is a proper scaling 
function. It takes very similar values in DP and MN, 
namely, z DP = 1.580745± 10" 6 while z MN = 1.53±0.07. 

We have performed numerical simulations of the SSW 
model for different q values starting from an initially 
pinned (p(0) = 1) interface. The exponents and 8 
have been measured by studying large system sizes (from 
L = 2 17 to L = 2 20 ) in order to minimize finite-size cor- 
rections, and by averaging over a small number of dif- 
ferent realizations (« 10) to further decrease statistical 
fluctuations. Conversely, it is sufficient to consider much 
smaller sizes (2 5 < L < 2 10 ) for a reliable estimate of 
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FIG. 2: Phase diagram of the SSW model. The depirming 
transition takes place along the solid line; the dashed line 
is the result of the approximate analytic mapping discussed 
in the text. In the shaded area, stationary pinned interfaces 
exist even though v w is smaller than the free interface velocity 
(=1/2). 



z that has been determined by looking for the optimal 
collapse of the various p(t, L) curves (in this case, how- 
ever, it has been necessary to average over a much larger 
ensemble of realizations (> 10 4 — 10 5 ). 

A complete summary of the resulting values of the crit- 
ical exponents is reported in Tabled Altogether, we find 
that for q < q*, the attracting force not only does leave 
the transition point unaffected, but also the universal- 
ity class of the critical behaviour remains of MN type 
[l8j . For q > q* , the transition line veers down, while the 
critical exponents signal a transition of DP type. 
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0.50(1) 
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0.276(5) 
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0.42975(5) 


0.16(1) 


0.27(1) 


1.5(1) 
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0.34895(5) 


0.17(1) 


0.276(5) 


1.5(1) 



TABLE I: Critical exponents of the SSW model. In paren- 
theses we report the estimated uncertainty on the last figure. 

The DP critical line can be best understood by ana- 
lyzing the SSW model in the vicinity of the point q = 1 
and v w — 0. Here, the dynamics is dominated by two 
slow mechanisms: (i) detachment of pinned sites during 
the asynchronous part of the rule; (ii) shrinking of the 
unpinned islands at the wall move. In comparison, the 
dynamics of detached regions between consecutive wall 
moves rapidly leads them to assume a perfectly triangu- 
lar shape with a maximal slope equal to ±1. Therefore, 
such islands correspond to the dead phase in DP, since 
their shape prevents the occurrence of any pinning in 
their interior. It is now convenient to divide between at- 
tached [hi = h w ) and detached {hi > h w ) sites, denoting 
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FIG. 3: The interface before (panel a) and after (b) the wall 
move, with the same settings as in Fig.Q The corresponding 
configurations read as ADADDDDDADDDDDDDADAD 
and DADADDDADADDDDDADADA (see text for the 
definition of A and D.) 

them with A and D, respectively. Correspondingly, the 
dynamics reduces to a simple probabilistic cellular au- 
tomaton: L/v w sites are first randomly selected, trans- 
forming each A into a D with probability 1 — q; next, the 
wall move amounts to transforming all A's into D's and 
all D's neighbouring an A into a D (see Fig. |3J. Apart 
from the peculiar presence of both an asynchronous and 
a synchronous part, this rule clearly belongs to the class 
of contact processes with an absorbing state (D) and, 
as such, it is expected to exhibit a DP transition pfl l2l| . 
The only relevant parameter is the ratio between the rate 
(1 — q) and the wall velocity v w , which corresponds to 
the slope of the critical line at q = 1. In fact, numeri- 
cal simulations of the automaton yield a "critical ratio" 
a c = lim^o v%j/(l — q) = 2.866 . . ., in very good agree- 
ment with the slope determined from direct simulations 
of the SSW model. 

At a finite distance from q = 1 , the automaton does no 
longer describe exactly the SSW dynamics, since fluctua- 
tions of the interface within an unpinned island can both 
induce a faster shrinking of the island and the generation 
of pinned sites in its interior. The first effect amounts to 
increasing by some factor a the average number of D's 
turned into A's at the island borders, so that the critical 
point would be determined by the equation 



In the configuration plotted in Fig. 1, the island bor- 
der identified by the letter C shifts by three units thus 
implying a = 3. The average value of a can be deter- 
mined by noticing that, at criticality, in the SSW model, 
the typical slope s of the interface inside an island must 
be such that its velocity coincides with the wall veloc- 
ity v^. If we then assume that the profile is a biased 
random walk with probabilities p u and 1 — p u of up and 
down moves, respectively, one finds, using simple com- 
binatorial considerations, a = 1 + 2(1 — p u )/p u . Since 
p u is simply related to the slope, s, by p u = (s + l)/2 
one finally obtains a = 1 + 2(1 — s)/(l + s). In- 
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serting this expression into © and eliminating s with 
the help of the relation = (1 — s 2 )/2, one obtains 
v c w {q) = 2y/l + 2a c (l - q) - a c (l - q) - 2. Although ap- 
proximate, this formula reproduces very accurately the 
DP critical line not only in the vicinity of q — 1 but also 
up to the multicritical point, where it touches the MN 
critical line (see Fig.|2Jl, providing a good approximation 
for its position as well (q max = 1 — 3/(2a c ) = 0.477. . .). 
Evidently, the quality of the theoretical formula implies 
that even close to the multicritical point, the sudden ap- 
pearance of pinned sites inside unpinned islands does not 
significantly modify the transition point. A deeper un- 
derstanding of this point is left to future investigations. 

The correspondence with MN and DP critical phenom- 
ena unveiled for small and large q values, respectively, 
is not sufficient to make predictions about the scaling 
behaviour in the vicinity of the multicritical point [l9|. 
There, in the absence of a convincing field-theoretic ap- 
proach, a chance for understanding how the two out-of- 
equilibrium critical phenomena may be connected to one 
another is offered by numerical investigation. However, 
even this is not a straightforward task, since three levels 
of criticality mix together: criticality of the free rough 
interface, criticality of the depinning transition and, fi- 
nally, that one connected with the MN-DP transition. 
While approaching q* from the left along the critical line 
v w = 1/2: t ne power-law decay of the density of pinned 
sites p{t) turns out to be first governed by the exponent 
8* « 1/2, which crosses over to Smn- As the crossover 
time appears to diverge when q — > q*, it can be safely 
stated that 8* characterizes the critical behaviour at q*. 



The same scaling is expected for q approaching q* from 
the right along the critical line, with a crossover from 8* 
to 8 dp- However, the difficulty of locating the critical 
line with a sufficient accuracy prevents an effective nu- 
merical verification. The other two critical exponents are 
[3* « 3/4 and z* w 3/2 at q* , so that both the [3* and 
the 8* values are intermediate between the corresponding 
MN and DP critical exponents, while z* is compatible 
with both zmn and zp>p. 

In conclusion, with reference to a simple microscopic 
model, we have shown that MN and DP can be differ- 
ent facets of the same wetting process. The connection 
between these two different universality classes strongly 
hints at the possibility that both out-of-equilibrium tran- 
sitions may be described within a single field theoretic 
approach. Some progress has been recently made in this 
direction in where the authors have found DP be- 
havior in a KPZ equation with attractive wall. 
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